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I. INTRODUCTION 

It is now clear that the origin of structure in the Uni- 
verse is a primordial density perturbation, existing al- 
ready when cosmological scales start to enter the hori- 
zon. Observation is consistent with the hypothesis that 
the density perturbation is perfectly adiabatic, Gaussian 
and scale-independent, but significant departures from 
this state of affairs is still allowed by the data. In partic- 
ular it is not excluded that the adiabatic density pertur- 
bation may be accompanied by a significant isocurvature 
density perturbation . 

Inflation provides a natural origin for the perturbation, 
since it converts the vacuum fluctuation of each light free 
scalar field into a classical scale-independent perturba- 
tion. One or more of these field perturbations may cause 
the primordial density perturbation. 

It is usually assumed that inflation involves a slowly- 
rolling field, dubbed the inflaton, whose value determines 
the end of inflation. The perturbation in the inflaton field 
cannot cause an isocurvature perturbation, but does in- 
evitably cause at some level an adiabatic perturbation. 
The usual hypothesis H is that the inflaton is solely re- 
sponsible for the observed adiabatic density perturba- 
tion. Under this 'inflaton hypothesis' significant non- 
Gaussianity is excluded in the usual one-field models 
{§,|). In multi-field models, where there is a family of 
possible inflaton trajectories curved in field space, sig- 
nificant non-Gaussianity is possible |^],f^] but apparently 
only at the expense of extreme fine-tuning of the ini- 
tial condition that specifies the inflaton trajectory. Any 
isocurvature density perturbation must be caused by the 
perturbation of some non-inflaton field. Under the infla- 
ton hypothesis this means that the isocurvature density 
perturbation (if present) depends on different physical 
parameters from the adiabatic density perturbation. As 
a result, an isocurvature perturbation of observable mag- 
nitude would require fine-tuning of the physical param- 
eters, or else some as-yet unforseen connection between 
them. 

An alternative hypothesis is that the adi- 

abatic density perturbation originates from the pertur- 



bation in some 'curvaton' field different from the infla- 
ton. In this scenario the adiabatic density perturbation 
is generated only after inflation, from an initial condi- 
tion which corresponds to a purely isocurvature pertur- 
bation p^j . □ The object of the present paper is to ex- 
plore the nature of the primordial density perturbation 
under this hypothesis. In the curvaton scenario, signifi- 
cant non-Gaussianity may easily be present because the 
curvaton density is proportional to the square of the cur- 
vaton field. Also, the curvaton density perturbation can 
lead, after curvaton decay, to isocurvature perturbations 
in the densities of the various components of the cosmic 
fluid. These, which we term 'residual' isocurvature com- 
ponents, are either fully correlated or fully anti-correlated 
with the adiabatic density perturbation, with a calcula- 
ble and generally significant relative magnitude. 

The paper is organized as follows. We deal in Section 
II with the adiabatic perturbation and its possible non- 
Gaussianity. In Section III we formulate the description 
of isocurvature perturbations, in a way which will allow 
us to analyse CDM, baryon and neutrino perturbations in 
a unified manner. In Section IV we calculate the residual 
isocurvature perturbations of cold dark matter (CDM) 
and baryons. In Section V we give a general formal- 
ism for describing the primordial neutrino isocurvature 
perturbation, taking into account for the first time the 
crucial issue of lepton number. Then wc use it to calcu- 
late the residual isocurvature neutrino perturbation. Our 
conclusions are summarised in Section VI. 



"This conversion mechanism has been considered also in the 
pre big bang scenario [HJli). In this scenario though, the 
required scale-invariant curvaton field perturbations will be 
generated only if the curvaton has anon-trivial coupling and 
for particular initial conditions pa,16|. 
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II. THE CURVATURE PERTURBATION 

A. The curvature perturbation and the primordial 
density perturbations 

From the viewpoint of observation the 'primordial' 
epoch is the one a few Hubble times before the small- 
est cosmological scale approaches the horizon. Taking 
that scale to enclose say 10 6 solar masses, the primor- 
dial epoch corresponds to temperature of order lOkeV 
which is after nucleosynthesis. Leaving aside the pos- 
sibility of a particle decaying after nucleosynthesis, the 
content of the Universe at the primordial epoch is there- 
fore rather well known. There are photons, practically 
massless neutrinos, baryons and (assuming it is already 
in existence) cold dark matter, with the radiation domi- 
nating the energy density. The corresponding primordial 
density perturbation has an adiabatic mode 
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which leaves the local ratio of number densities unper- 
turbed. Non-adiabatic or isocurvature modes are defined 
by 
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In this section we focus on the adiabatic mode, which is 
known to be the dominant one responsible for structure 
formation, and return to the possible isocurvature modes 
in the next section. 

The primordial adiabatic density perturbation is asso- 
ciated with a spatial curvature perturbation. Following 
[ p7| , ^8[ we define the curvature perturbation ( on spatial 
slices of uniform density p with the line element 



d£ 2 = a 2 (l + 2C)<5ydiz;M2 J 
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The quantity C, is related to the density perturbation, Sp, 
and curvature perturbation, ip 7 on a generic slicing (using 
the sign, convention of Rcf. [^9|) by the gauge- invariant 
formulalO 



C = -ip - HAt 
P 



(6) 
(7) 



^In general a gauge corresponds to a definite slicing and 
threading of spacetime, but in this paper only the former 
is relevant so that 'gauge- invariant' can be taken to mean 
'independent of the slicing'. 



where At is the displacement of the generic slicing from 
the uniform-density slicing. On super-horizon scales £ 
is practically identical with the curvature 1Z defined on 
slices orthogonal to comoving worldlincs. The quan- 
tity £ w 1Z is useful on such scales because it is time- 
independent po| , ^l| , p2|Jl^ j3]l provided that the pressure 
perturbation is adiabatic, meaning that SP = c 2 8p, 
where the adiabatic sound speed c 2 = P/p. This is guar- 
anteed if there exists a universal equation of state P{p). 
If the pressure perturbation is not adiabatic, then £ « 1Z 
changes according to the equation pOp^,pj 
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where <5P na d = SP — c 2 Sp. 

In the conventional inflaton model for the origin of 
structure purely adiabatic perturbations are generated 
due to quantum fluctuations in the single scalar field driv- 
ing inflation. Thus the curvature perturbation, £*, cal- 
culated shortly after Hubble-exit (k = aH) determines 
the curvature perturbation until that scale re-enters the 
Hubble scale during the subsequent radiation or matter- 
dominated era. By contrast, we are interested in a 
scenario where the curvature perturbation generated on 
large-scales during inflation in the very early universe is 
negligible. 

In the rest of this section we describe the generation 
of the large-scale curvature perturbation, in the curva- 
ton scenario, amplifying the outline given in the original 
paper ]9| (see also [PU ). We deal with the simplest ver- 
sion of the curvaton scenario, where the curvature per- 
turbation is caused exclusively by a single 'curvaton' field 
distinct from the inflaton field. 



B. Generating the curvaton field perturbation 

During inflation the Hubble parameter H is assumed 
to be slowly varying, en = H/H 2 <C 1. The inflaton (if 
it exists) is supposed to produce a negligible curvature 
perturbation. Demanding that it is (say) less than 1% of 

the observed value implies || V* < 2 x I0 15 GeV, which 
in turn implies [^) that the primordial gravitational waves 
will have no detectable effect on the CMB anisotropy. 
Conversely, the detection of an effect would rule out the 
curvaton hypothesis (an "anti-smoking gun" ) |2^,|2^] . 

The curvaton field a is supposed to be practically free 
during inflation, with small effective mass ( | V^-o- 1 <C H 2 
where a subscript a denotes d/da). It follows that on 
super-horizon scales there is a Gaussian perturbation 
with an approximately scale- independent spectrum given 
by 



(9) 



where the * denotes the epoch of horizon exit k = aH. 
(The normalization of the spectrum || is such that Pg a 
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specifies the typical magnitude of a spatial fluctuation in 
5a on a physical scale a/k.) The spectral index specifying 
the slight scale-dependence is given by |p|,p5t 

n - 1 = d ln(^v)/d In k = 2r, aa - 2e H , (10) 

where = V arT /3H 2 . 

After the smallest cosmological scale leaves the hori- 
zon, the curvature perturbation remains negligible until 
after the curvaton starts to oscillate. As a result, the cur- 
vaton field evolves in unperturbed spacetime. To follow 
the evolution it is assumed that the curvaton field has 
no significant coupling to other fields or, to be more pre- 
cise, that the effect of any coupling can be integrated out 
to give a possibly time-dependent effective potential V. 
With these assumptions, the curvaton field smoothed on 
the smallest cosmological scale evolves along each comov- 
ing worldline according to the unperturbed field equation 



a + 3Ha + Vv = . 
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Making the first-order approximation (5(V CT (x, i)) « 
V aa (t)Sa(x,t), its inhomogeneous perturbation satisfies 
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(12) 



Assuming only that the evolution of 5a is linear, the 
Gaussianity and (slight) scale-dependence of the origi- 
nal quantity are preserved. The fractional perturbation 
5a(x)/(7 does not evolve at all if the potential is either 
sufficiently flat (so that a and 5a are both constant) or 
quadratic (so that a and 5a satisfy the same equation). 

The field remains over-damped until the Hubble pa- 
rameter falls below the curvaton mass m a . The curvaton 
field will then start to oscillate about its vacuum value 
(taken to be a — 0) with an amplitude which decreases 
with time. Even if the potential of the curvaton field is 
not quadratic, after a few Hubble times we can make the 
approximation V ~ ^m 2 a 2 . The fractional field pertur- 
bation then has some constant value, 
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where the factor q is time-independent because the oscil- 
lation is around a quadratic minimum of the potential. 
In particular, if the effective potential for the curvaton is 
quadratic (or sufficiently flat) throughout its evolution, 
q = l. Using Eq. ([)[), the spectrum of the fractional field 
perturbation at this stage is 



The energy density in the oscillating field is 
p CT (x,£) = m 2 a a 2 (x,t) , 
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where ct(x, t) is the amplitude of the oscillation. The 
perturbation in p a depends on the curvaton field per- 
turbation through both a linear and a quadratic term. 



Assuming, for the moment, that the linear term domi- 
nates (which we shall see is demanded by the data) we 
have 



Pa 



6a 
2 — 
a 



,' 5a 
2q[ — 
a 



(16) 



We shall return to consider the possible contribution from 
the quadratic term in section II D . 



C. Generating the curvature perturbation 

So far we have reached the epoch just after the Hub- 
ble parameter falls below the curvaton mass, and the 
curvaton oscillation starts. At this stage, it is supposed 
that the dominant energy density comes from radiation. 
The curvaton, though, is supposed to be fairly long-lived, 
while decaying comfortably before nucleosynthesis. So 
long as the decay-rate is negligible (r <C H) 1 we have 
p a oc a~ 3 and p r oc a~ 4 , leading to p a /p r oc a. It is this 
increase in the relative curvaton energy density which 
generates the curvature perturbation. 

To analyse the generation of the curvature perturba- 
tion (and also, later on, to discuss possible isocurvature 
perturbations produced by the curvaton decay) , it is con- 
venient to consider the curvature perturbations Q corre- 
sponding to the separate components of the energy den- 
sity. These are defined on slices of uniform pi, corre- 
sponding to the gauge-invariant definition 
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In particular, evaluating £ CT for the curvaton on unper- 
turbed (ip — 0) hypersurfaces when the curvaton starts 
to oscillate, we have 
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The total curvature perturbation (]?]) can then be writ- 
ten as the weighted sum 



C = (1 - /)Cr + fC 



(19) 



where the relative contribution of the curvaton to the 
total curvature is given by 
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Thus the curvaton perturbation, C, a , can initially be de- 
scribed as an isocurvature perturbation since / — > 0, i.e., 
Pal Pr ~ * 0, in the early-time limit. 

Until the effect of curvaton decay becomes significant, 
the curvaton and radiation densities each satisfy their 
own energy conservation equation pi = —3H(pt + Pi). 
In this regime, each Q is constant on super-horizon 
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scales Jig] . The evolution of C on these scales is due 
solely to the change of / in Eq. ( |l9| ) , which yields 



C = /(Cr-Cr) = #/(!-/) (Ca-Cr 



(21) 



It is intuitively easy to see how ( changes with time as the 
curvaton density grows relative to the radiation, / > 0. 
It is straightforward to check that Eq. (^lj) is consistent 
with Eq. (||) where the non-adiabatic pressure perturba- 
tion is given by 
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An adiabatic perturbation Sp r /p r = Sp a / 'pa corresponds 
to the special case 



C = Cr = Co 
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The curvaton scenario corresponds to the case where 
the curvature perturbation in the radiation produced at 
the end of inflation is negligible, Cr ~ 0. In the ap- 
proximation of sudden decay |J , Cr and Co- both remain 
constant up until decay and the curvature perturbation 
at decay is therefore 



C ~ /dec Co 
~ nJdccq 



(24) 
(25) 



where /d GC is / at the decay epoch, conventionally defined 
in terms of the decay rate by H^ cc = T. We assume that 
after decay, the pre-existing radiation is cither insignif- 
icant, or else thermalizes with the decay products with 
the possible exception of CDM. Since we are assuming 
that there is only one curvaton field, this is sufficient to 
ensure that C on cosmological scales will retain the same 
value until the primordial epoch. 

Going beyond the sudden decay approximation, we de- 
fine a number r by 
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where C is evaluated well after the epoch of curvaton de- 
cay and Co- is evaluated well before this epoch. In the 
limit where the curvaton completely dominates the en- 
ergy density before it decays, r = 1. In other words, 
the sudden-decay approximation becomes exact in this 
limit. The reason is that in this case the curvaton and 
its decay products constitute a single fluid with a definite 
P(p), so that they have constant curvature perturbation 
which is equal to Co-- If the curvaton does not dominate, 
one has to resort to numerical calculation of the coupled 
perturbation equations [Eg], for which one finds 
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The prediction of the curvaton model for the spectrum 
of the curvature perturbation is 



n = 3^/0/0 



(30) 



with Vg a i a given by Eq. (|l4|). The COBE measurement 
of the CMB quadrupole anisotropy requires 



VI (COBE) = 4.8 x 10" 
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If the curvaton dominates the energy density before it 
decays (r = 1) this implies the following amplitude for 
the perturbations of the curvaton field 
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In other words, the typical field (and density) pertur- 
bation on cosmological scales must be of order 10 -4 in 
that case. If the curvaton does not dominate, we need at 
least r > 10~ 4 to get a spectrum of the correct mag- 
nitude (since the typical density perturbation can be 
at most of order 1). We shall now see how the obser- 
vational bound on non-Gaussianity actually requires a 
much higher value. 



Possible non-Gaussianity of the curvature 
perturbation 



D. 



From Eqs. ( [15|) and (|27j), the curvature perturbation 
depends on the curvaton field perturbation through both 
linear and a quadratic terms. So far we evaluated only 
the linear contribution, which gives a Gaussian contri- 
bution to the curvature perturbation. However, if the 
quadratic term in the density perturbation of the curva- 
ton field is not negligible, the total curvature perturba- 
tion will have a non-gaussian {x 2 ) component. 

The relative magnitude of the quadratic term is con- 
ventionally specified || by a number /nl (NL mean- 
ing 'non-linear'), which nominally determines the non- 
Gaussian contribution to the Bardeen potential accord- 
ing to the formula 
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The connection with the Bardeen potential is only nomi- 
nal, because the relation between $ and C is taken || to 
be 
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This relation is actually valid only on super-horizon scales 
after matter domination, the correct relation at the pri- 
mordial epoch being more complicated and involving the 
relative neutrino density |||. 
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Using Eq. ( |l5| ) we have 
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Thus, using Eqs. (|33| ) and (|34j) , the prediction of the 
curvaton hypothesis is 



(36) 



At this point we have to remember that first-order cos- 
mological perturbation theory is being assumed. Second- 
order metric perturbations will also give a correction 
$( 2 )(x), presumably with typical magnitude 
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It follows that the estimate Eq. (pq ) of /nl will be valid 
only in the regime /nl 3> 1 (and that smaller values of 
/nl cannot even be defined, unless $( 2 )(x) is itself the 
square of a Gaussian quantity). In other words, the va- 
lidity of the estimate Eq. ( |36| ) requires that the curvaton 
contributes only a small fraction of the energy density 
just before it decays. In the opposite case that the cur- 
vaton dominates the density before it decays, the non- 
Gaussianity calculated according to linear theory is lost 
in the "noise" of the unknown second-order correction to 
cosmological perturbation theory. 

Now we compare the prediction with present and fu- 
ture data on the CMB anisotropy. A recent analysis of 
the COBE data §7| yields |/ NL | < 2 x 10 3 , implying the 

constraint r > 6 x 10~ 4 . From Eqs. ( |30|) and ( p"ll) we see 

i 

that this bound is equivalent to V^ a , g < 0.1. Starting 
in 2003, data from the MAP satellite p^l will either de- 
tect non-Gaussianity or give |/nl| < 2C l|| corresponding 
to r > 0.06. If non-Gaussianity is detected by MAP it 
will clearly be above the noise of the second-order cor- 
rection. Looking further ahead, the PLANCK satellite 
will either detect non-Gaussianity or give |/nl| < 5 ||] 
corresponding to r > 0.2. Apparent non-Gaussianity at 
the bottom end of this range would however have to be 
checked against the second-order order correction to $ 
to make sure that the correction is negligible. 

We have seen that the curvaton hypothesis can easily 
give significant non-Gaussianity. The reason is that the 
predicted curvature perturbation is proportional to the 
perturbation of the curvaton density, a quantity which 
depends on the square of the curvaton field. The per- 
turbation in the curvaton field is supposed to be Gaus- 
sian because it is supposed to have negligible interaction, 
making its Fourier components uncorrelated which is the 
definition of Gaussianity. But the corresponding pertur- 
bation in the curvaton density is a linear combination 
of the curvaton field perturbation and its square. From 
a theoretical viewpoint the square could even dominate 
p9|,p[, though as we have seen this is ruled out by the 
data. 



This is in sharp contrast with the situation for the in- 
flaton hypothesis, in which the curvature perturbation is 
purely linear in the inflaton field perturbation. The in- 
flaton can usually be treated as a practically free field, 
making its perturbations and the curvature practically 
Gaussian. In particular, for the usual single-field mod- 
els, where there is an essentially unique inflaton trajec- 
tory, the self-interaction of the inflaton field is the only 
relevant one and is kept small by the flatness condi- 
tions e,)) € 1. As a result single- field models give Q] 
/nl| = |2e — 2rj\ < 0.1 which is lost in the noise from 
the second-order correction to cosmological perturbation 
theory.El In two-field models, where there is a family of 
possible inflaton trajectories curved in field space, the in- 
flaton field perturbation at the end of inflation will be a 
linear combination IBOtl 
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where 4>* is the inflaton field at the time of horizon exit 
and x* is t ne field orthogonal to it. By definition the 
slope of the potential vanishes in the x direction but the 
higher derivatives might in principle be present corre- 
sponding to a self-interaction big enough to generate sig- 
nificant non-Gaussianity. In particular there might be a 
cubic term Mx 3 in the potential |H| (with the unper- 
turbed x* — 0) or something more complicated This 
setup might generate significant non-Gaussianity, but it 
obviously requires a special choice of the inflaton trajec- 
tory. Even when such a choice is made initially, it might 
be destabilized by the quantum fluctuation before cos- 
mological scales start to leave the horizon. 

We conclude that a detection of non-Gaussianity by 
MAP would strongly suggest that the primordial curva- 
ture perturbation is produced by a curvaton field. Such 
a detection would imply that the density of the curvaton 
before it decays reaches no more than 6% of the total. In 
the opposite case that the curvaton dominates before it 
decays, it gives a practically Gaussian curvature pertur- 
bation, just as in the inflaton scenario. 



III. ISOCURVATURE PERTURBATIONS 

Thus far we have concentrated upon how the curva- 
ton introduces a large-scale curvature perturbation. Now 
we consider the isocurvature perturbations that may be 
imprinted in different particle species after the curvaton 
decay. 

We adopt the 'separate universe' viewpoint, implicit 
in practically all discussions of perturbations on super- 
horizon scales. At each epoch, it is assumed that each 
comoving region with size much bigger than the Hubble 



'The bound 0.1 comes from the observational bound on the 
spectral index in the inflaton scenario, n = 2rj — 6e. 
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distance looks locally like some unperturbed (Robertson- 
Walker) universe. If these 'separate universes' are all 
identical, the cosmological perturbations are said to be 
adiabatic. Perturbations in matter fields then vanish on 
slices of uniform energy density, and the curvature per- 
turbation C is the only thing that needs to be specified to 
determine the evolution of the perturbations after hori- 
zon entry. 

If the 'separate universes' are not identical there are 
isocurvature perturbations, so-called because they evolve 
independently of the curvature perturbation on large 
scales which can therefore be taken to vanish when con- 
sidering them. One way of specifying a generic isocur- 
vature perturbation, 5x, is to give its value on uniform- 
density slices, related to its value on a different slicing by 
the gauge-invariant equation 



rii'4) 

X J Sp=0 
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(39) 



For a set of fluids with energy density pi, the isocurvature 
perturbations are instead conventionally defined by the 
gauge-invariant quantities 



-3H 
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Pi 
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In terms of the individual curvature perturbations de- 
fined in Eq. ( p"7j ) this becomes 



(41) 



For fluids which are 'decoupled' (in the sense that there 
is no energy transfer) and which have a definite equa- 
tion of state Pj(/9j), the Q's are constant on superhori- 
zon scales [[18| and hence so are the isocurvature per- 
turbations (JH|). In this language, the curvaton density 
perturbation before it decays corresponds to a constant 
isocurvature perturbation S ar = 3 (Co- — Cr) which in the 
curvaton scenario (Cr = 0) reduces to S ar = 3C CT . 

At the 'primordial' epoch, before the smallest cosmo- 
logical scale approaches the horizon, there are four 'de- 
coupled' fluids namely the CDM, baryons, photons and 
neutrinos (taken to be massless). The conventional def- 
initions Eqs. (§)-(§) of the three isocurvature perturba- 
tions correspond to 



Sb — 3 (Cb - Cy) 

^cdm = 3 (Ccdm Cy) 
Sy = 3 (Ciy — C7 ) ) 



(42) 
(43) 
(44) 



In the standard single-inflaton scenario it is impossible 
to introduce isocurvature perturbations on large scales at 
the primordial epoch from the purely adiabatic pertur- 
bations existing after the end of inflation |l8| . But in the 
curvaton scenario the non-adiabatic nature of the curva- 
ton perturbation after inflation {S ar ^ 0) means that it 
is possible for the curvaton to leave isocurvature pertur- 
bations after the curvaton decays even on super-horizon 



scales, which we term 'residual' isocurvature perturba- 
tions. 

In addition to the energy densities we shall need to con- 
sider 'number' densities n c dm, n-B, til and nn. They are 
defined, respectively, as (i) the number density of CDM 
particles, (ii) the density of baryon number B, (iii) the 
density of lepton number L and (iv) the densities of the 
three individual lepton numbers Li (i = e, p or r). By 
the 'primordial' era any significant lepton number will be 
carried entirely by neutrinos, making nn the difference 
between neutrino and anti-neutrino number densities. In 
the unperturbed Universe each number density is pro- 
portional to a~ 3 so long as the corresponding quantity is 
conserved in a comoving volume. 

For each of these number densities, it will be useful 
to consider the curvature perturbations Ci on slices of 
uniform n, , 



-1> 
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In a homogeneous universe the conserved number density 
rii obeys the evolution equation 



3Hrn = . 



(47) 



Allowing for large-scale perturbations about the strictly 
homogeneous background, we obtain the local evolution 
equation for the first-order perturbation 



5n,i + 3HSr 



3riiip = , 



(48) 



assuming all spatial gradients (e.g., divergence of the par- 
ticles velocity field) are negligible on large scales. Com- 
bining equations (|45|-|4^) we obtain the evolution equa- 
tion for the curvature perturbation 



Ci = o. 



(49) 



Hence we see that Q defined in Eq. ( j45| ) remains con- 
stant on large-scales so long as the particle number is 
conserved. 

Equations ( f45| ) and (|9|) are equivalent to the state- 
ment that the perturbation in the fractional number den- 
sity (Srii/rii) is conserved on flat (ip = 0) spatial hyper- 
surfaces. The perturbations in the number densities are 
best defined on the flat slicing, because on this slicing 
the expansion rate with respect to coordinate time is un- 
perturbed jlSfl , leading to the constancy of the fractional 
perturbations. 

This is a significant extension to the case of interacting 
fluids of the result found in Ref. Jjl| for the constancy 
of Ci for non-interacting fluids on large-scales. It will 
prove a powerful tool to describe the possible genera- 
tion of isocurvature perturbations after curvaton decay 
in what follows. 
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IV. RESIDUAL ISOCURVATURE MATTER 
PERTURBATIONS 

A. Residual isocurvature CDM perturbation 

Without making any assumption about the nature of 
the CDM, we take the epoch of its creation as the one 
after which its particle number n c d m is conserved. There 
are several candidates for the CDM particle, such as the 
axion, the lightest supersymmetric particle (LSP), a sta- 
ble massive particle (wimpzilla) or primordial black holes. 
For axions the epoch of creation corresponds to the tem- 
perature ~ 1 GeV at which the effective axion mass rises 
above the Hubble parameter, and for the conventional 
LSP it corresponds to the freeze-out temperature roughly 
of order 10 GeV. For very massive particles (wimpzillas) 
p2| , |33[ or black holes |34|,[55| the epoch of creation may 
be the end of slow-roll inflation corresponding to an en- 
ergy scale up to 10 16 GeV. In the case of wimpzillas 
it may instead [p6| be the epoch of thermal inflation cor- 
responding to perhaps ~ 10 GeV. The CDM might 
also be created as an out-of-equilibrium decay product of 
the inflaton, the curvaton or some other particle. 

After CDM creation, conservation of the CDM par- 
ticle number ensures that Ccdm defined by Eq. (|||) is 
constant on super-horizon scales . When the CDM be- 
comes non-relativistic with constant mass, Ccdm reduces 
to Ccdm defined by Eq. ([l7j). In most of the examples 
mentioned above, this occurs at the epoch of creation. It 
could happen though that that CDM created from out-of- 
cquilibrium decay is initially relativistic. Also, the axion 
mass (induced by the QCD instanton) increases after the 
epoch of creation, becoming constant only when the tem- 
perature falls to 100 MeV. In any case, Ccdm will certainly 
have reduced to Ccdm by the primordial epoch. 

To evaluate the primordial isocurvature perturbation 
Scdm, we assume that the primordial C7 — C correspond- 
ing to a zero or small neutrino isoc urva ture perturbation 
Cr/ — Cy> postponing until Section V A the possibility of 
a significant neutrino isocurvature perturbation. We can 
then re-write Eq. (E3I) as 
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(50) 



As we have mentioned, several different mechanism 
have been considered for creating the CDM. In most 
cases, the creation mechanism does not involve any quan- 
tity with an isocurvature perturbation, which means that 
at creation the CDM has no isocurvature perturbation. 
To put it differently, the 'separate universes' are in these 
cases identical at the epoch of CDM creation, which 
means that the number density n c d m at creation depends 
only on the local energy density. 

The most usual way of obtaining CDM which has an 
isocurvature perturbation at creation is to suppose that 
it consists of axions, whose production in each region of 



space corresponds to the onset of the axion field oscilla- 
tion in that region. If the axion field has an isocurvature 
perturbation, a CDM isocurvature perturbation will then 
be produced at the time of creation. However, the mag- 
nitude of the isocurvature CDM perturbation produced 
in this way is unrelated to that of the curvature pertur- 
bation. In other words, there is no reason why CDM 
originating from an oscillation of the axion field should 
have a significant isocurvature density at the time of its 
creation. The same is true for all of the other CDM pro- 
duction mechanisms that have been considered so far. 

In what follows we will assume that the CDM num- 
ber density at the epoch of creation depends only on the 
local energy density. We expect this to be valid in the 
absence of any non-adiabatic pressure perturbation, i.e., 
where the local density also determines the local pres- 
sure. The only exception will be the case that the CDM 
is produced directly by the curvaton decay, which will 
require separate treatment. 

In the curvaton scenario where the CDM is created af- 
ter the curvaton has decayed (or in the inflaton scenario) 
the CDM isocurvature perturbation will be zero at the 
'primordial' epoch. The basic reason is that the 'sepa- 
rate Universes' in these cases are identical. To proceed 
more formally, the assumption that there is initially no 
isocurvature density perturbation means that the CDM 
number density at creation is uniform (5n c dm — 0) on 
slices of uniform density (where ip = — £)• Hence we may 
evaluate the gauge-invariant expression (^5|) on a surface 
of uniform density to obtain Ccdm = C with both sides 
constant. Going forward to the primordial epoch, when 
the CDM energy is conserved so that Ccdm = Ccdm we 
find that the primordial CDM isocurvature perturbation, 
•Scdm given by Eq. (^), vanishes. 

The situation is dramatically different if the CDM is 
created before the curvaton decays, or if it is created 
by the curvaton decay itself. In these cases, the process 
of curvaton decay creates a 'residual' isocurvature CDM 
perturbation whose properties are closely tied to the cur- 
vature perturbation. 

Consider the case that the CDM is created well before 
the curvaton decays and well before the curvaton comes 
to dominate the energy density (/ C 1 in Eq. (p0|)). 
In this regime the Universe is practically unperturbed 
and therefore Ccdm is practically zero. At the primordial 
epoch this gives the residual CDM isocurvature pertur- 
bation in Eq. (|3|) as 



5 c dn 



-3C. 



(51) 



In the notation of Ref. g] this corresponds to a maximum 
correlation between curvature and CDM-isocurvature 
perturbations. The power spectra of the two perturba- 
tions have the same spectral index and the isocurvature 
perturbation has an amplitude three times larger than 
the adiabatic one. Such a large correlated perturbation 
is ruled out by current observations which require P]37|] 
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< 1.5 at 95%c.l. 



(52) 



Our conclusion is that in the curvaton scenario, CDM 
cannot be created before the curvaton decays and while 
the total curvature perturbation remains small, ( <C rC, a . 
In particular CDM creation just after inflation ends, such 
as might occur in the case of wimpzillas or black holes, 
is incompatible with the curvaton scenario. 

If the CDM is created before the curvaton decays, but 
when the curvaton density has become non-negligible, 
/ ~ 1 in Eq. (^) , there will be a significant non-adiabatic 
pressure perturbation, <5P na d given by Eq. (p2|), at the 
epoch of creation and we can no longer assume that the 
initial CDM number density is unperturbed on uniform 
density slices. Instead the actual number density at cre- 
ation will depend on the mechanism by which the CDM 
is created, and we leave a detailed investigation of the 
different cases for future work. 

Finally, we consider the case that the CDM is created 
by the decay of the curvaton itself. The epoch of CDM 
creation then corresponds to the epoch when the curva- 
ton decay is complete. The resulting local CDM density 
is then a fixed multiple of the curvaton number density 
well before decay. The fractional perturbations are thus 
equal and hence 



Ccdm Ccr 



(53) 



Using the definition of r in Eq. (26), this gives at the 
primordial epoch 



Scdm — ^ [ ) C • 



(54) 



This is negligible if the curvaton comes to dominate be- 
fore it decays (r = 1), the physical reason being that 
the curvaton perturbation becomes an adiabatic curva- 
ture perturbation, so it cannot leave behind any residual 
isocurvature perturbation. But if the curvaton decays 
before it dominates, r«l, there will be large isocurva- 
ture perturbations in the CDM, incompatible with exist- 
ing observational constraints. Note that in the notation 
of Ref. the curvature and isocurvature perturbations 
are anti-correlated, in which case the observational limit 
becomes 



< 0.2 at 95%c.l. 



(55) 



much stronger than that for correlated isocurvature per- 
turbations. Observational limits on the amplitude of the 
anti-correlated isocurvature perturbations thus require 
the curvaton to effectively dominate, r > 0.9, if the CDM 
is created by the curvaton decay. 



B. Residual baryon isocurvature perturbation 

Baryon number B may be created directly, in which 
case we take the epoch of baryogenesis to be the one after 



which B is conserved. Alternatively B may be created 
through the conversion of L at the electroweak transition. 
In the latter case B — L is conserved, and for the present 
purpose we may take the epoch of baryogenesis as the 
one when B — L is created. Depending on what type of 
mechanism operates, the epoch of baryon creation in this 
sense may be anywhere from the end of inflation to the 
electroweak transition. To keep the notation simple we 
take the relevant quantity to be B from now on, with the 
understanding that B is to be replaced by B — L if that 
is the relevant quantity. 

Until the QCD transition at T ~ 100 MeV, B rep- 
resents the difference between the abundances of the 
typically relativistic particles and anti-particles carrying 
baryon number, and cannot usefully be associated with 
an energy density. In this situation we need to use £b the 
curvature perturbation on slices of uniform ub (given 
in Eq. (pi|)). This quantity is constant on large scales 
at all times after baryon creation, and it reduces to £g 
(Eq. (|l7j)) after the QCD transition when baryon number 
is carried by non-relativistic nucleons and nuclei. 

The theoretical situation for the baryon isocurvature 
perturbation is similar to the one for CDM. We make the 
reasonable assumption that the baryon number density 
at the epoch of baryogenesis depends only on the local 
energy density, unless baryogenesis comes the curvaton 
decay itself. If baryogenesis occurs after the curvaton 
decays, there will be no primordial baryon isocurvature 
perturbation. If instead baryon number is generated well 
before curvaton decay and well before the curvaton dom- 
inates, there will be a large residual baryon isocurvature 
perturbation, 



s B = -K- 



(56) 



In contrast with the CDM case, such a perturbation is 
still marginally compatible with current observational 
data The effect on the CMB angular power spec- 
trum of a baryon-isocurvature perturbation is essentially 
the same as that of a CDM-isocurvature perturbation but 
the size of the effect is diminished by a factor Vis/ ^ c dm 
due to the smaller density of baryonic matter |39| , so the 
constraint in Eq. (52) becomes 



Sb 
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< 1.5 



at 95%c.l. 



(57) 



An isocurvature perturbation of the form given in 
Eq. ( [56] ) will either be ruled out or observed in the near 
future, providing in the latter case a smoking gun for the 
curvaton scenario. 

The final possibility is that baryon number is produced 
by the out-of-equilibrium decay of the curvaton itself, as 
is the case if we identify the curvaton with the right- 
handed sneutrino of ref. ]3§| ] . Then we have at the pri- 
mordial epoch, in the sudden decay approximation, 
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(58) 
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Current observational limits on the amplitude of such 
an anti-correlated baryon-isocurvature mode require that 
the curvaton must dominate the density at decay r > 0.6, 
if the curvaton itself is responsible for baryogenesis. 



V. NEUTRINO ISOCURVATURE 
PERTURBATIONS 

A. The general formalism 

In this section we discuss the possibility of a neutrino 
isocurvature density perturbation. As far as we are aware 
our treatment is the first one taking into account the 
crucial issue of lepton number. 

The era of thermal equilibrium for neutrinos ends just 
before nucleosynthesis, which almost certainly means 
that this era determines the abundances of neutrinos and 
anti-neutrinos at the later 'primordial' epoch. (We are 
not in this paper considering the possibility that the cur- 
vaton decays after nucleosynthesis, and for the present 
purpose we discount too the possibility that it decays 
between neutrino decoupling and nucleosynthesis.) Tak- 
ing that for granted, there is no neutrino isocurvature if 
the lepton number at decoupling is negligible, because 
the primordial abundance of the neutrinos is then deter- 
mined entirely by the photon temperature |3!| . 

We therefore consider the case of non-zero lepton num- 
ber density TiLi, with i =e, [i or r. While the neutrinos 
are effectively massless and in equilibrium, with temper- 
ature T v this corresponds to occupation numbers 



fi(E) = [exp(E/T v t &) + !]- 



(59) 



for neutrinos and anti-neutrinos with energy E. This 
gives the following expressions for the total energy density 
Pi of neutrinos and anti-neutrinos, and the lepton number 
density nn equal to the difference between the number 
densities of neutrino and anti- neutrinos [El: 
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(60) 
(61) 

(62) 
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(64) 
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In the usual case that the lepton asymmetry is negligi- 
ble, these expressions start to become valid at reheating, 
and hold with T v = T 7 until positron annihilation, after 
which 



T„ = (4/ll)*T 7 



(66) 



In the case of significant lepton asymmetry that we are 
considering, there is significant neutrino mixing and the 
expressions become valid only when T v — T 7 falls to a 
few MeV. The subsequent evolution is more complicated 
than in the usual case so that the thermal distribution 
Eq. (|S9|) is not precisely maintained pQ] , but following 
the usual practice we make the approximation that it 
is maintained so that Eqs. ([39])-(|67J) are all valid after 
positron annihilation. From now on we take the expres- 
sions to refer to that era. The asymmetry parameters £j 
are then constant since the neutrinos are travelling freely 
leading to n j oc oc 1 /a? . 

Big-bang nucleosynthesis (BBN) and large-scale struc- 
ture (LSS) constrain £ e and |£ p | 2 + |£ T | 2 - If the favoured 
large mixing angle (LMA) solution to the solar neutrino 
problem is correct, neutrino oscillations ensure that the £j 
have a common value £ 0] . In that case, the BBN/LSS 
constraint is J4§] -0.01 < £ < 0.07. With this con- 
straint, unperturbed lepton number densities are almost 
certainly too small to have any observable effect on the 
CMB anisotropy or large-scale structure (LSS). We shall 
see, though, that this need not be the case for the residual 
isocurvature perturbation. 

Final confirmation of the LMA is expected in a 
few months from the Kamland reactor experiment E3j , 
but we shall nevertheless allow independent asymme- 
try parameters in our analysis. In that case BBN and 
CMB/LSS together give constraints on = |£ T | (as- 
sumed equal for convenience) and £ e which are much 
weaker, namely |£ M , r | < 2.5 and |£ e | < 0.30 @. Note 
that these limits were obtained assuming a purely adia- 
batic primordial perturbation spectrum and hence serve 
only as a rough guide in the case of a correlated curvature 
and neutrino isocurvature perturbations. 

For each individual neutrino species we can calcu- 
late the curvature Q on uniform-/^ hypersurfaces from 
Eqs. @ and (|6T 
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and the curvature Q on uniform-nii hypersurfaces given 
from Eqs. (45) and (O) as 
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(68) 



where A^ = (60/7)5, and = 1 + 3(&/tt) 2 . Only when 
8£i = do these two quantities coincide. 

As we have noted, a neutrino isocurvature perturba- 
tion due to the first term in Eqs. (^) and (|8|), repre- 
senting a perturbation in the neutrino to photon tem- 
perature ratio, is extremely unlikely since it would have 
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to be generated during the extremely short era between 
neutrino decoupling and nucleosynthesis. We here take 
the first term to be zero. 

Following the path that we trod for the CDM and 
baryon number perturbations, we define the isocurvature 
perturbations in the number densities of the three lepton 
numbers, 



Sni 



5n~. 



= 3(5 - Cr) • 



(69) 
(70) 



They determine the perturbations in the asymmetry pa- 
rameters, 



(71) 



where the final equality is valid in the regime (Ci/V) 2 <C 
1. These expressions refer to the primordial era, and 
are valid in the early Universe only back to the epoch of 
positron annihilation at T ~ MeV. While they are valid, 
the lepton numbers Li are conserved and the curvature 
perturbation Q on slices of uniform nn is constant. At 
early times though, corresponding to T bigger than a few 
MeV, neutrino mixing becomes significant and the indi- 
vidual lepton numbers are not defined. Instead there is 
only the total lepton number density and its perturba- 
tion, 



Sl 
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n Li 



5n L 
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(72) 
(73) 
(74) 



The total lepton number L is well-defined in the early 
Universe after the epoch of lepton number creation, 
which we take to be the one after which L is conserved. 
While L is conserved, the curvature perturbation £l on 
slices of uniform til is constant. 

The three primordial quantities m and their pertur- 
bations Si are not in general determined by and Sl. 
They are however determined if the LMA solution to the 
solar neutrino problem is correct because the three quan- 
tities are then equal 



1 

nLi = g«L 
Si = Sl 



(75) 
(76) 



If one or more neutrino masses are big enough to give a 
significant amount of dark matter, it will be necessary to 
insert the 6£i into the initial occupation number Eq. (59), 
propagating this initial condition forward with evolution 
equations. In that case the initial condition requires a 
knowledge of the three individual isocurvature perturba- 
tions Si. Here we make instead the opposite approxi- 
mation of massless neutrinos. Then, the only required 



initial condition is the isocurvature perturbation in the 
total neutrino energy density which is specified by the 
derived quantity S v , defined by Eq. (0). This quantity 
may be calculated using 



taking pi from Eq. (|60j) and Q from Eq. (|6 
is 



(77) 



The result 
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(78) 
(79) 



where the final approximation is valid for (£i/7r) 2 <C 1. 

As far as we are aware, this is the first time that ex- 
pressions for the neutrino isocurvature perturbation have 
been given in the most realistic case where it is deter- 
mined by the (perturbed) lepton asymmetry. The ef- 
fect of the neutrino isocurvature perturbations has never 
been tested against observational data while including 
non-zero chemical potential. 



B. Induced matter isocurvature perturbations 

Now we address an issue that is relevant for calcu- 
lating the amplitude of CDM and baryon isocurvature 
perturbations in the presence of a neutrino isocurvature 
perturbation. In the presence of the latter, the equality 
C7 — Cv = C breaks down, and instead we have 

C = (1 - R,)(- ( + R„( v = (- ( + ^S u , (80) 

where R v = p v /{p v + /%) is the fraction of the final ra- 
diation density in neutrinos, and R v w 0.41 for ^ < 1. 

The amplitude of the CDM and baryon isocurvature 
perturbations, calculated assuming £ — £ 7 , acquire an 
additional term 

"Scdm/B = 3(Ccdm/B — C7) = 3(Ccdm/B ~ C) + RyS v . 

(81) 

Hence for CDM/baryons created after curvaton decay (or 
in the inflaton scenario), for which Ccdm/s — Cj we now 
have a non-zero isocurvature matter perturbation 



S c 

dm/B — Ri/Su , 



(82) 



for CDM/baryons created before curvaton decay, with 
Ccdm/s = 0, we have 



S c dxa/B — ~ 3£ + R V S V 



(83) 



while for CDM/baryons created by curvaton decay, with 
Ccdm/s = Cr, we have 
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<S cdm / B = 3 (—^ J C + R V S V . (84) 

If S u <C C, then £ 7 can be identified to high accuracy with 
C, and the changes in 5 c d m /B induced by the neutrino 
isocurvature perturbations represent a small correction. 

These changes do not correspond to any change in the 
evolution of the CDM or baryons in the early Universe. 
They would be avoided if we worked instead with the 
quantities 5 c dm/B = 3(C c dm/B — C)> that define the per- 
turbations on slices of uniform total density (Eq. (p9[)). 
These quantities are independent of S„ so that, for exam- 
ple, they always vanish if the CDM or matter is created 
after curvaton decay with a density depending only on 
the local energy density. In a similar way, the quan- 
tities Si = 3(£i — C) reflect the early Universe situation 
more directly than the quantities Si , as we shall see when 
discussing leptogenesis from curvaton decay. In all cases 
though, the unhatted quantities are those commonly used 
at the 'primordial' starting point for the forward 

evolution of the perturbations. 

C. The residual neutrino isocurvature perturbation 

Now we wish to discuss the residual neutrino isocur- 
vature perturbation in the curvaton scenario, along the 
same lines as for the CDM and baryons. To do this, we 
assume that the LMA solution is correct so that the pri- 
mordial lepton number densities have a common value, 
determined by the total lepton number density that is 
conserved in the early Universe. Since the lepton asym- 
metry is small in that case, the neutrino isocurvature 
perturbation is given to high accuracy by Eq. ( |79| ) , 

S*«y(|) (85) 

where the common asymmetry parameter £ satisfies the 
nucleosynthesis constraint |£| < 0.07. 

To evaluate Sl we can follow closely the previous dis- 
cussions of CDM and baryon number. We take the epoch 
of creation of lepton number L to be the one at which 
this quantity starts to be conserved. We assume that the 
lepton number isocurvature perturbation at creation is 
zero or negligible, except in the case that lepton num- 
ber is created by the curvaton decay. If lepton number 
is created after curvaton decay, Sl = and there is no 
neutrino isocurvature perturbation. 

If lepton number is created well before curvaton decay 
and well before the curvaton dominates the density, 

Sl = -3C 7 , (86) 

giving 

S„ = -^(£)V (87) 



For |£| < 0.07 and using Eq. (j8(i|), we have 

|5„|<0.01C. (88) 

This is very small and may never by observable. 

Finally, if lepton number is created by out-of- 
equilibrium curvaton decayB 

^-sf^V (89) 

If the curvaton does not dominate before it decays, this 
gives 

*-?(H3(§)'c 

\s v \ < o.oi i-Ic<i5C, (91) 

where the final inequality comes from the current 
bounds on the non-Gaussianity parameter, /nl, given 
in Eq. (^). This upper bound represents a large anti- 
correlated neutrino isocurvature component which may 
already be incompatible with current data, though this 
possibility is yet to be tested against observations. If 
MAP further constrains the non-Gaussian parameter, 
/nl, (rather than detecting non-Gaussianity) the bound 
will become |S„| < 0.2£. 

We note that the neutrino isocurvature velocity mode 
considered by Bucher et al |39| is hard to generate in 
a curvaton scenario, since we would expect any velocity 
perturbation left on large scales after the curvaton decay 
to be suppressed by factors of order (k/aH) 2 with respect 
to the (almost scale-invariant) density perturbation, i.e., 
generated only by spatial gradients in the curvaton field. 

VI. CONCLUSIONS 

In this paper we have reviewed the mechanism by 
which a curvaton field can generate a large-scale den- 
sity perturbation after inflation, and gone on to inves- 
tigate the nature of the primordial perturbation that is 
produced. In many cosmological models the curvaton 
scenario can reproduce the purely adiabatic and Gaus- 
sian density perturbation (curvature perturbation) that 
is familiar from the inflaton scenario. In another regime 
though, it can give large non-Gaussianity, and/or a large 
'residual' isocurvature perturbation which is completely 
correlated with the curvature perturbation and can in 
principle be present in any or all of the baryonic matter, 
the CDM and the three neutrino species. 



' The exact formula is Sl = 3f(l— r)/r where Sl = 3(£l — C), 
but the difference is negligible in the regime (£/tt) 2 <C 1. 
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The non-Gaussianity of the curvature perturbation £ 
arises if the curvaton fails to dominate the energy density 
before it decays. The perturbation is described by a x 2 - 
distribution, whose non-Gaussianity is parameterised by 



/nl — 



4r 



(92) 



where r is approximately given by the curvaton density 
just before decay, as a fraction of the total. In 2003, 
results from MAP satellite will either detect this non- 
Gaussianity or show that the curvaton density is at least 
a few percent of the total when it decays. 

The residual isocurvature perturbation in CDM or 
baryonic matter arises if the CDM or baryon number is 
created either before curvaton decay, or by the curvaton 
decay itself. If it is created significantly before curva- 
ton decay and before the curvaton dominates the energy 
density, the isocurvature perturbation is given by 



"Scdm — 3£ , 



(93) 



and similarly for baryons. This large, correlated isocur- 
vature perturbation is ruled out by observation for CDM. 
In other words, CDM cannot be created significantly be- 
fore curvaton decay and before the curvaton dominates. 
This is a strong constraint on the cosmology, implying for 
example that the curvaton scenario is inconsistent with 
the CDM creation just after the end of inflation. For 
baryons the perturbation amplitude is close to current 
observational limits and will be ruled out or detected in 
the near future, providing in the latter case a smoking 
gun for the curvaton scenario. 

If the CDM or baryon number is created by the curva- 
ton decay itself, 
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(94) 



and similarly for baryons. Unless r is close to 1, this large 
anti-correlated contribution is ruled out by observation 
for both CDM and baryons. In other words, the curvaton 
decay cannot create CDM or baryons unless the curvaton 
dominates the energy density before it decays. 

For the neutrino isocurvature perturbation we have 
presented an analysis which, for the first time, includes 
the crucial issue of lepton number, and is relevant for 
any cosmology. Although general formulae are given, we 
focus on the case where the ratio of the neutrino and 
photon temperature is unperturbed, because the late de- 
coupling of neutrinos makes it very hard to see how it 
could be otherwise. We also focus on the case that the 
e, \x and r lepton number densities are equal, which will 
be ensured by mixing in the early Universe at least if the 
large mixing angle solution to the solar neutrino problem 
is correct. 

With these assumptions, the isocurvature energy den- 
sity perturbation S v of massless neutrinos is related to 
the lepton number density isocurvature perturbation Sl 
by 



S - = T 1 

7 \TT 



(95) 



where £ is the lepton asymmetry parameter, related to 
the lepton number density by £ — 4.02(?il/?i 7 ). If the 
lepton number density is of order the baryon number 
density as is usually supposed, |£| ~ 10~ 9 and the neu- 
trino isocurvature perturbation will be completely unde- 
tectable. Observationally though, the lepton density is 
subject only to the nucleosynthesis bound |£| < 0.07, 
which may allow a detectable perturbation. However 
neutrino isocurvature perturbations in the presence of 
a significant lepton asymmetry have never been tested 
against observations. 

In the curvaton scenario, the formulas for the resid- 
ual isocurvature perturbation Sl are the same as those 
for the CDM and baryon perturbations. If leptogene- 
sis occurs after curvaton decay there is no residual neu- 
trino isocurvature perturbation. If it is significantly be- 
fore curvaton decay and before the curvaton dominates, 
Sl = —3( and |S„| < 0.01C which may never be ob- 
servable. The interesting result comes in the third case, 
that the curvaton decay itself causes leptogenesis. Then 
Sl = 3((1— r)/r)Cj and if the curvaton does not dominate 
before decay we have a neutrino isocurvature perturba- 
tion whose magnitude is related to the non-Gaussianity 
parameter /nl = 5/4r, 
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(96) 



If the present bound on non-Gaussianity, |/nl| < 2000, 
is saturated this permits a huge effect, presumably al- 
ready ruled out by observation, and even the expected 
MAP bound on |/nl| < 20 will still allow a significant 
effect. Observational bounds on /nl and <5>„ in this sce- 
nario should be obtained jointly, taking into account the 
correlation of these two quantities with each other and 
with the curvature perturbation 

In summary, we have shown that different curvaton 
scenarios offer a number of distinctive observational pre- 
dictions which may be tested by forthcoming experi- 
ments. 

Note added: A related paper by Moroi and Taka- 
hashi |Tl]| discussing the case of matter isocurvature per- 
turbations in the curvaton scenario appeared while this 
work was in progress. 



ACKNOWLEDGMENTS 

We are grateful to Kostas Dimopoulos, Luca 
Amendola and Martin Buchcr for useful discus- 
sions. This work was supported in part by PPARC 
grants PPA/G/S/1999/00138, PPA/G/S/2000/00115 
and PGA/G/O/2000/00466, and by EU grant HPRN- 
CT-2000-00152. DW is supported by the Royal Society. 



12 



[1] R|, Trotta, A, Kiazuclo and K. Durrcr, F hys. Key. Lett. 87 



[28] http://map.gsfc.nasa.gov/ 
[29] D. H. Lyth, Phys. Rev. D 45 (1992) 3394. 
[30] C. Gordon, D. Wands, B. A. B assett and R. Maartens 
Phys. Rev. D 63, 023506 (2001) ferXiv:astro-ph/000913ll 
[31] N. Bartolo, S. Matarrese and A. Riotto, u:Xiv:hep- 



2 31301 (2001) arXiv:astro-ph/01CB 0T7[ 



ph/0112261 



[2] L. Amendola, C. Gordon, D. Wands and M. Sasaki, P hys. 

Rev. Lett. 88, 211302 (2002) |arXiv:astro-ph/0107088| ]. 
[3] A. R. Liddle and D. H. Lyth, "Cosmological inflation and 

large-scale structure," Cambridge, UK: Univ. Pr. (2000) 

400 p. 

[4] A. Gangui, F. Lucchin, S. Mat arrese and S. Mollerach , 
Astrophys. J. 430, 447 (1994) fcrXiv:astro-ph/9312033l; 
A. Gangui Phys. Rev. D 50, 3684 (1994) ]arXiv:astr& 
ph/9406014|; L- M. Wang and M. Kamionkowski, Phys 



Rev. D 61, 063504 (2 000) [|arXiv:astro-ph/9907431 
A. Gangui and J. Martin, arXiv:astro- ph/9908009| . 



E. Komatsu and D. N. Spergel, Phys. Rev. D 63 (2001) 
063002. 

D. S. Salopek, P hys. Rev. D 45 (1992) 1 139; F. Bernardeau 

and J. P. Uzan, |arXiv:hep-ph/0207295 . 

N. Bartolo, S. Matarrese and A. Riott o, Phys. Rev. D 65 

(2002) 103505 |arXiv:hep-ph/0112261 |. 

A. D. Linde and V. Mukhano v, Phys. Rev. D 56, 535 

(1997) |arXiv:astro-ph/9610219| ]. 

D. H. Lyth and D. Wa nds, Phys. Lett. B 524, 5 (2002) 
[ arXiv:hep-ph/011000S ]. 

T. Moroi and T. Takah ashi, Phys. Lett. B 522, 215 (2001) 
| arXiv:hep-ph/0110096 ]. 



T. Moroi and T. Takahashi, |arXiv:hep-ph/0206026 . 
S. Mollerach, Phys. Rev. D 42, 313 (1990). 
K. Enqvist and M. S. S loth, Nucl. Phys. B 626, 395 (2002) 
|arXiv:hep-ph/0109214]. 



V. Bozza, M. Gasperin i, M. Giovannini and G. Veneziano, 



arXiv:hep-ph/0206131 



E. J. Copeland, R. Easther and D. Wands, Phys. Rev. D 



56, 874 (1997) [a.rXiv:hep-th/970108: 



J. E. Lidsey, D. Wands and E. J. Cope land, Phys. Rept. 
337, 343 (2000) JarXiv:hep-th/9909061 |. 
J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. 
Rev. D 28, 679 (1983). 

D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle , Phys. 
Rev. D 62, 043527 (2000) ]arXiv:astro-ph/0003278|[. 



V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, 

Phys. Rept. 215, 203 (1992). 

J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

D. H. Lyth, Phys. Rev. D 31, 1792 (1985). 

J. Martin and D. J. Sc hwarz, Phys. Rev. D 57, 3302 (1998) 

| arXiv:gr-qc/9704049 |. 

J. Garc ia-Bellido and D. Wand s, Phys. Rev. D 53, 5437 
(1996) |arXiv:astro-ph/951102s|l. 
D. H. Lyth, |arXiv:hcp-ph/0205266 . 

D. Wands, N. Bartolo, S. Matarrese and A. Riotto, 



arXiv:astro-ph/0205253 



K. A. Malik, D. Wands and C. Ungarelli, in preparation. 
E. Komatsu, B. D. Wandelt, D. N. Spergel, A. J. Ban- 
day and K. M. Gorski, Astrophys. J. 566 (2002) 19 



arXiv:astro-ph/0107605 


. M. G. Santos et ai, Phys. Rev. 


Lett. 88 (2002) 241302 


arXiv:astro-ph/010758£]. 



[32] D. H. Lyth, D. Roberts fc M. Sm ith, Phys. Rev. D 57 



7120 (1998) | |arXiv:hep-ph/960944l| ] 
[33] D 



J. Chung, E W. Kolb and A. Riott o, Phys. Rev. D 59, 
023501 (1999) |arXiv:hep-ph/9802238 1 . 
L. Randall, M. Soljacic and A. H. Guth, Nucl. Phys. B 



472, 377 (1996) |arXiv:hep-ph/951243S| ]. 
J. Garcia-Bellido, A. D. Linde and D. Wan ds, Phys. Rev. 
D 54, 6040 (1996) |arXiv:astro-ph/9605094 1 . 
L. Hui and E. D. Stewart, Phys. Rev. D 60, 023518 
(1999) jirXiv:liep-ph/9812345ll; G. 
B 459, 489 (1999) |arXiv:hep-ph/9905204| : 
G il, E. J. Copeland, J . Gray, A. Lukas and M. Plumacher, 



Dyali, Phys. Lett. 
M. Bastero- 



[34] 
[35] 
[36] 

[37] L. Amendola, private communication. 

[38] K. Hamaguchi, H. Mu rayama and T. Yanagid a, Phys. Rev. 

D 65, 043512 (2002) ferXiv:hep-ph/010903C |. 
[39] M. Bucher, K. Moodley and N. Turok Phys. Rev. D 62, 

083508 (2000) |arXiv:astro-ph/990423l| 
[40] A. D. Dolgov' kXiv:hcp-ph/0202T22; ; G. G. Raffelt, 
|arXiv:astro-ph/020722c| ; K. N. Abazajian, J. F. Bea- 
com and N. F. Bell, Phys. Rev. D 66 (2002) 013008 

Phys. Rev. D 66 



arXiv:hep-th/020104C 



JarXiv:astro-ph/0203442|; Y. Y . Won; 



(2002) 025015 |arXTyThep-ph/0203180 |. 
[41] A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, 
G. G Raffelt and D. V. Sem ikoz, Nucl. Phys. B 632 (2002) 
363 ]arXiv:hep-ph/0201287 l. 
[42] M. Orito, T. Kajino, G. J Mathews and Y. Wang , Phys. 

Rev. D 65, 123504 (2002) ]arXiv:astro-ph/0203352| . 
[43] Talk at Neutrino 2002, for transparencies 



http://neutrino2002.ph.tum.de/ 



13 



